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ABSTRACT 

NGC 1851 is an intriguing Galactic globular cluster, with multiple stellar evolution¬ 
ary sequences, light and heavy element abundance variations and indications of a 
surrounding stellar halo. We present the first results of a spectroscopic study of red 
giant stars within and outside of the tidal radius of this cluster. Our results identify 
nine probable new cluster members (inside the tidal radius) with heliocentric radial 
velocities consistent with that of NGC 1851. We also identify, based on their radial 
velocities, four probable extratidal cluster halo stars at distances up to ~3.1 times the 
tidal radius, which are supportive of previous findings that NGC 1851 is surrounded 
by an extended stellar halo. Proper motions were available for 12 of these 13 stars 
and all are consistent with that of NGC 1851. Apart from the cluster members and 
cluster halo stars, our observed radial velocity distribution agrees with the expected 
distribution from a Besangon disk/N-body stellar halo Milky Way model generated 
by the Galaxia code, suggesting that no other structures at different radial velocities 
are present in our field. The metallicities of these stars are estimated using equiva¬ 
lent width measurements of the near infrared calcium triplet absorption lines and are 
found, within the limitations of this method, to be consistent with that of NGC 1851. 
In addition we recover 110 red giant cluster members from previous studies based on 
their radial velocities and identify three stars with unusually high radial velocities. 

Key words: globular clusters: general; globular clusters: individual: NGC 1851; stars: 
kinematics and dynamics; techniques: spectroscopic; techniques: radial velocities 


1 INTRODUCTION 

Until relatively recently globular clusters (GCs) were con¬ 
sidered as examplars of single stellar populations (SSPs), 
with all the cluster stars having the same age, helium abun¬ 
dance and overall [Fe/H] metallicity. In reality all GCs that 
have been thoroughly studied have multiple stellar popula¬ 
tions (with the possible exception of Ruprecht 106, but the 
sample size was smallVillanova et al. 2013). Evidence comes 
from many photometric studies showing multiple branches 
on colour-magnitude diagrams (CMDs), and high resolution 
spectroscopic studies showing that all globular clusters ex¬ 
hibit significant star-to-star abundance variations in light 
elements (such as C, N, O and Na). In addition a few clus¬ 
ters (typically the most massive) also show overall [Fe/H] 
metallicity variations and variations in neutron capture ele¬ 
ments. 

Globular clusters are, therefore, more complex then 
originally thought, and in fact contain multiple generations 
of stars (for a review see Gratton et al. 2012). This was first 


discovered for the more massive globular clusters and for 
the brightest red giant branch (RGB) and asymptotic giant 
branch (AGB) sequences. Now many globular clusters with 
a wide range of masses have been shown to exhibit these fea¬ 
tures, and various studies have shown multiple sequences in 
the fainter main sequences (MS), sub giant branches (SGB) 
and horizontal branches (HB). 

The abundances of light elements of stars in GCs ex¬ 
hibit various relations, in particular C-N and O-Na anticor¬ 
relations. The origin of these patterns is believed to be hot 
hydrogen burning in either AGB stars, fast rotating massive 
stars, massive binaries or supermassive stars (Gratton et al. 
2012). However there are indications that the C-N and O- 
Na anti-correlations are not perfectly coupled. Pancino et al. 
(2010) measured CN and CH bandstrengths for MS stars 
in 12 GCs with well-characterized RGB abundances and 
found that the frequency of C-poor, N-ric.h MS stars does not 
match the frequency of O-poor, Na-rich RGB stars reported 
for those same clusters in Carretta et al. (2009). They sug- 
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gested that enrichment in C and N may not match enrich¬ 
ment in O and Na because the CN cycle, the full CNO bi¬ 
cycle and the Ne-Na cycle are each responsible for a different 
part of the light element abundance pattern, and each re¬ 
quires a different temperature to operate. 

NGC 1851 is an unusual Galactic globular cluster. It is 
one of four globular clusters (along with NGC 1904, NGC 
2298 and NGC 2808) all physically located within a sphere 
of radius 6 kpc (Bcllazzini et al. 2004). It has been proposed 
that they form a group associated with the putative Canis 
Major dwarf galaxy stream (Martin et al. 2004). NGC 1851 
exhibits the common features of multiple evolutionary se¬ 
quences and star-to-star abundance dispersions in the light 
elements, but it also has some features not ordinarily found 
in other globular clusters, as detailed below. 


1.1 Photometry 

Stetson (1981) detected peculiarities in the HB with a dense 
clump of stars at the red end and a hint of another clump 
at the blue end. Walker (1992) confirmed the split nature 
of the HB in the cluster, clearly shown in the right panel 
of Fig. 8. A double SGB was announced in Milone et al. 
(2008b); they postulated two populations with similar over¬ 
all [Fe/H] metallicity but an age difference of ~1 Gyr, while 
Cassisi et al. (2008) suggested as an alternative two coeval 
populations with different C+N+O abundance mixtures. 
Han et al. (2009) also detected a split RGB using UI pho¬ 
tometry. These studies strongly suggest that a second gener¬ 
ation of higher overall [Fe/H] metallicity stars formed after 
an early episode of metal enrichment. 


1.2 Spatial extent 

There have been a number of photometric and spectroscopic 
studies (detailed below) that indicate that NGC 1851 has a 
surrounding stellar halo outside the tidal radius. The nature 
of this halo is not yet clear - they may be stars stripped by 
tidal or disk shocking. Alternatively, modelling has shown 
that it is possible that NGC 1851 is the nuclear star cluster 
of a now disrupted dwarf galaxy, with the halo stars repre¬ 
senting the remnant of that galaxy. 

The tidal radius rt calculated by different dynamical 
models such as King (King 1966), Wilson (Wilson 1975) 
and power-law templates can vary considerably, and various 
studies have adopted different values. The current version 
of the Harris 1996 (2010 edition) catalogue does not tab¬ 
ulate rt values, but instead lists the structural parameters 
core radius r c = 0.09 arcmin and central concentration c = 
1.86 arcmin from McLaughlin & van der Marel (2005). The 
central concentration c is given by: 

c = lo 3y (!) 

From this the tidal radius rt is calculated as 6.5 ar¬ 
cmin. The Wilson (1975) model rt as calculated 
by McLaughlin & van der Marel (2005) is 44.7 arcmin. 
Carballo-Bello et al. (2012) looked at the number density 
profiles of 19 Galactic globular clusters and found power-law 
templates to be a better representation for about 2/3 of their 
sample, including NGC 1851 for which their calculated tidal 


radius is 11.7 arcmin. Allen et al. (2006) calculated the Ja¬ 
cobi (or Roche) radius as 17 arcmin, although this is strongly 
dependent on the adopted Galaxy gravitational potential 
and cluster orbital parameters. For results discussed in this 
introduction, we quote the tidal radius adopted by the au¬ 
thors in the paper where applicable. For our own study we 
adopt the value of Carballo-Bello et al. (2012) of 11.7 ar¬ 
cmin. 

Wide field photometric star-count analyses carried out 
by Leon et al. (2000) show a symmetric stellar structure 
out to ~20 arcmin and extensions to ~45 arcmin that 
are roughly aligned towards the Galactic center. Photo¬ 
metric studies by Olszewski et al. (2009) also show indi¬ 
cations of a stellar halo extending up to 75 arcmin, ^6.5 
times beyond their adopted tidal radius of 11.7 arcmin. This 
halo has a mass of ~1 per cent of the dynamical mass of 
NGC 1851, but there was no evidence for any substructure or 
tidal stream formation. Carballo-Bello & Martfnez-Delgado 
(2010) reported a distinct MS of a metal-poor stellar over¬ 
density with a sky projected size >15.7° in the area sur¬ 
rounding NGC 1851 and NGC 1904. 

A possible tidal tail was reported by Sollima et al. 
(2012) in a spectroscopic survey of 107 stars from 12 ar¬ 
cmin to 33 arcmin around the cluster centre. They not only 
found a significant fraction of stars outside their adopted 
tidal radius of 11.7 arcmin, but also an interesting hint of a 
kinematically cold stellar stream, albeit with a peak in the 
heliocentric radial velocity (V r ) distribution of ~180 kms^ 1 , 
compared with 320.5 ± 0.6 kms' 1 Harris 1996 (2010 edi¬ 
tion) for NGC 1851. Marino et al. (2014) carried out a spec¬ 
troscopic analysis of 23 stars in the halo of NGC 1851 and 
found stars at distances up to ~2.5 times their adopted tidal 
radius of 11.7 arcmin. However they did not find convincing 
evidence for the tidal tail of stars reported by Sollima et al. 
(2012) at ~180 kms -1 . Kunder et al. (2014) reported nine 
RAVE stars up to 10° away from NGC 1851 that may be 
associated with the cluster. 

Bekki & Yong (2012) simulated the dynamical evolu¬ 
tion of a dwarf galaxy with a stellar nucleus in the tidal 
field of the Galaxy. They showed that a diffuse stellar halo 
around NGC 1851 would be produced after stripping the 
dark matter halo and stellar envelope, and that the rem¬ 
nant nucleus was consistent with an object like NGC 1851. 
They also showed that mergers of globular clusters in host 
dwarf galaxies were possible, which in the case of NGC 1851 
could explain its multiple stellar populations. 

1.3 Chemistry 

NGC 1851 exhibits the typical GC light element star-to- 
star abundance dispersions and relations that are believed 
to result from multiple star formation episodes. Even for 
GC systems with multiple stellar populations it is unusual 
in being one of a handful of clusters that show evidence (de¬ 
tailed below) of variations in overall [Fe/H] metallicity and 
C+N+O abundances. It shows some similarities to cj Cen- 
tauri, which exhibits large star-to-star abundance variations 
in all elements (Lee et al. 1999 and Bedin et al. 2004), but 
also some similarities to clusters like NGC 6752 that have 
no [Fe/H] variation but show light element abundance vari¬ 
ations (Yong et al. 2013). However, some results have been 
conflicting, especially the question of whether there is a scat- 
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ter in the sum of the C+N+O abundances and whether this 
sum correlates with the metallicities or heavy-element abun¬ 
dances. 

Eight bright giant member stars analysed by 
Yong & Grundahl (2008) (hereafter YG08) showed 
abundance variations of the s-process elements Zr and 
La that were correlated with A1 and anticorrelated with O. 
There was a hint that these s-process element abundances 
clustered around two values, but the sample size was small. 
An analyis of high-resolution spectroscopic data for four of 
the YG08 RGB stars was reported in Yong et al. (2009). 
These stars showed a large C+N+O abundance variation 
that correlated with the light elements Na and Al, and also 
with the s-process elements Zr and La. 

The largest extant high resolution spectroscopic study 
of NGC 1851 stars, Carretta et al. (2010, 2011) (hereafter 
CIO), obtained spectra of 124 RGB stars and found a 
measurable [Fe/H] spread of 0.06 to 0.08 dex. A division 
into metal-rich and metal-poor groups based on [Fe/H] 
and [Ba/H] also shows many differences in various abun¬ 
dance correlations and anti-correlations. They proposed that 
NGC 1851 is the product of a merger of two globular clus¬ 
ters, with an age difference of ~1 Gyr and different Fe and 
a-capture element abundances, that formed in a now dis¬ 
rupted dwarf galaxy. 

There have been several other spectroscopic studies of 
NGC 1851 stars. Lardo et al. (2012) obtained C and N abun¬ 
dances in 64 SGB and MS stars and found significant anti¬ 
correlations in [C/H] and [C/N] and significant variations 
in C, N and total C+N content when divided into fainter 
and brighter samples. Villanova et al. (2010) studied a sam¬ 
ple of 15 red giant branch (RGB) stars in NGC 1851. Their 
analysis did not show variations in iron-peak elements or to¬ 
tal C+N+O content. In contrast, Yong et al. (2015) found 
a large (~0.6 dex) spread in the sum of the C+N+O abun¬ 
dances in AGB and RGB stars, providing support for the 
idea of coeval SGB populations with different C+N+O con¬ 
tent. 

While there has been a large-scale study of C and N in 
SGB and MS stars (Lardo et al. 2012), and a smaller study 
in RGB stars (Villanova et al. 2010), C and N abundances 
have not been measured for a large sample of RGB stars, 
and not for stars with measured light element abundances. 
This study was specifically aimed at expanding the set of 
abundance information available for the RGB stars in the 
CIO studies. In particular we aimed to carry out a thorough 
survey of C and N abundances for the same stars, which can 
be determined from moderate-resolution blue spectra. Com¬ 
bined with the proton-capture (O, Na, Mg, Al), a-capture 
(Si, Ca, Ti) and s-process (Y, Zr, Ba, La) element abun¬ 
dances obtained in their study, this will enable an assess¬ 
ment of whether the C-N and O-Na anticorrelations are tied 
tightly together in NGC 1851, or whether the relative con¬ 
tributions of the different proton-capture cycles to cluster 
self-enrichment can vary. In addition, we expect to settle 
the questions of whether there is a scatter in the sum of 
the C+N+O abundances and whether the sum correlates 
with the metallicities or heavy-element abundances deter¬ 
mined by CIO (see the conflicting results mentioned above 
from Carretta et al. 2011, Yong et al. 2009, Villanova et al. 
2010 and Yong et al. 2015). A second paper (in preparation) 


will address the determination of C and N abundances and 
discuss these questions. 

Because of the large field-of-view of AAOmega, we were 
also able to investigate the reports detailed above in Subsect. 
1.2 that there is an extended halo of stars around NGC 1851. 
This paper describes the observations and first data reduc¬ 
tion and analysis of the sample obtained as well as measure¬ 
ments of V r for all observed stars, and we assess the evidence 
for a stellar halo around NGC 1851. 


2 SAMPLE SELECTION, OBSERVATIONS 

AND INITIAL DATA REDUCTION 

We used the AAOmega spectrograph on the Anglo- 
Australian Telescope (AAT) at Siding Spring Observatory 
in two degree field (2dF) fibre positioner multi-object spec¬ 
trograph (MOS) mode to obtain the spectra 1 . The target 
list was made up of the 124 RGB stars from the CIO study 
and four of the giant stars analysed by YG08 (the main tar¬ 
gets of our study), plus a third field sample made up of stars 
with similar 2MASS Skrutskie et al. (2006) photometry to 
that of NGC 1851 RGB stars. The latter were selected with 
the aim of allowing for serendipitous discovery of new cluster 
members, extratidal cluster halo stars or of field stars with 
unusual abundance patterns. Observing details are listed in 
Table 1. 

The red channel spectra were used for determination 
of stellar heliocentric radial velocities lV r ) and metallicities 
[Fe/H]. We used the standard 2dF data reduction (2dFdr) 
pipeline, as well as the IR AF task splot to remove skylines. 

A total of 1876 spectra were obtained. Spectra with 
low signal-to-noise ratio (SNR) and/or bad columns or other 
artefacts were discarded leaving a total of 1730 spectra for 
the analysis. Because of the selection algorithm of the 2dF 
configuration software, 541 stars were observed more than 
once in different pointings. In total the useable spectra of 
1149 unique stars were obtained, including 106 stars from 
the CIO study and the four stars from YG08. In the vicinity 
of the CaT in the red channel, the median SNR of spectra 
obtained was ~28. Spectra of the star used as a template 
for calculation of V r and of one of our probable new member 
stars are shown in Fig. 1 to illustrate the quality of our data 
and the similarity of the spectra in our observations. 


3 HELIOCENTRIC RADIAL VELOCITIES 

We used the IRAF task fxcor to derive the V r of ob¬ 
served stars using the CIO star [CLG2011] 30965 (2MASS ID 
05142874-4003159) in our sample as a template star. Mean 
V r values were calculated for the 541 stars with multiple 
observations by weighting their fxcor V r values with their 
respective fxcor V r errors. 

Multiple observations of the same target were also used 
to estimate the V r errors for stars. An estimate of the stan¬ 
dard deviation of V r (cr v ) was calculated using the simpli¬ 
fied statistics for small number of observations methods of 
Dixon & Massey (1951). These a v for multiple observations 

1 Manuals and technical details at 

http/www.aao.gov.au/2dF/aaomega 
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Table 1. Observing details. 


Date 

Start Time 
(UTC) 

Field name 

Magnitude 

limits 

Exposures 
(number x s) 

Seeing 

(arcsec) 

Airmass 
(mid exposure) 

2012 Dec 17 

11:30 

bright-fl 

V < 14.5 

3 x 540 

1.4 

1.10 

2012 Dec 17 

12:30 

faint-fl-long 

V > 14.5 

2 x 2700 

1.4 

1.02 

2012 Dec 17 

14:21 

faint-fl-short 

V > 14.5 

1 x 2700 

1.3 

1.05 

2012 Dec 17 

15:44 

faint-f2-long 

V > 14.5 

2 x 2700 

1.5 

1.11 

2012 Dec 18 

10:51 

faint-f2-short 

V > 14.5 

1 x 2700 

1.6 

1.14 

2012 Dec 18 

15:34 

faint-f3-long 

V > 14.5 

2 x 2700 

1.4 

1.24 



8300 8400 8500 8600 8700 8800 

Wavelength (A) 


335 


330 


325 

'uj 

| 320 
>" 

O 315 


310 


305 


305 310 315 320 325 330 335 

C10V r (km s"* 1 ) 



Figure 1. Example spectra. Top panel: red channel spectrum 
of the star used as the template for fxcor V r measurements - 
CIO Simbad reference [CLG2011] 30965; 2MASS identification 
05142874-4003159. Bottom panel: red channel spectrum of one of 
our probable new member stars - 2MASS identification 05141398- 
4000194. The Calcium II triplet absorption lines are indicated. 
Note: flux scales are different on each plot. 

were then compiled as a function of the median continuum 
level of the spectrum between the two strongest CaT lines. 
The values of the mean a v for each bin for stars with multi¬ 
ple observations were then used as the errors AV r for stars 
with single observations as a function of their continuum 
level. For stars with N observations, AV r was calculated by 
dividing a v by -\/N. 

We compared our V r values for the stars in common 
with the CIO study. Fig. 2 shows the CIO V r of all 106 stars 
recovered in our study compared to our V r measurements 
of the same stars. Six of the 106 stars (5.6 per cent) had 
measured V r more than ±lcr different from those listed in 
CIO. Binarity is a possible explanation for these differences 
- binary stars will have variable V r , but will likely still have 
single-line spectra since our targets are red giant stars. The 
fraction of stars with V r more than one ila different is con¬ 
sistent with the estimated binary fraction of stars in glob¬ 
ular clusters (2 per cent for the most massive clusters to 
20 per cent for the less massive ones, with a mass weighted 


Figure 2. V r of all observed CIO stars compared to our V, mea¬ 
surements of the same stars. The solid black line is the linear 
best fit to our data, the solid green line shows the one-to-one re¬ 
lationship and the dashed green lines show ±1 <t of the one-to-one 
relationship. 


average of ~5 per cent, Milone et al. 2008a). It is also worth 
noting that three of these stars are among the fainter stars 
in our sample, so measurement of their V r may be affected 
by the lower SNR of their spectra rather than by binarity. 

Selection as a probable new member was based on the 
V r range of the recovered CIO and YG08 groups. The min¬ 
imum and maximum measured V r for stars in these groups 
was 308.0 kms -1 and 328.9 kms^ 1 respectively. We found 
13 stars in our field sample with V r within this range. 

Fig. 3 shows the V r distribution of observed stars. The 
left panel shows all the observed stars, with the predicted 
V r distribution of the Milky Way model generated by the 
Galaxia code (Sharma et al. (2011), see Sect. 6) overplot¬ 
ted as a grey histogram. The distribution shows two promi¬ 
nent peaks: i) a peak with a mean V r ~40 kms -1 con¬ 
sisting of 1015 stars with a spread in radial velocity from 
-100 kms -1 < V r < 200 kms -1 that we identify as Galac¬ 
tic disc or halo field stars (non cluster members), and ii) 
a peak with a mean V r ~320kms~ 1 containing 123 stars 
that we identify as members of NGC 1851. In the follow- 
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Figure 3. Left panel: V r distribution of all observed stars with the predicted V r distribution of the Milky Way model generated by the 
Galaxia code (Sharma et al. 2011) overplotted as a grey histogram. Right panel: Expanded V r distribution of cluster member and cluster 
halo stars. Blue denotes known CIO members, green known YG08 members and red are stars that are probable new cluster members or 
cluster halo stars. The vertical dashed line shows the Harris 1996 (2010 edition) catalogue V r of 320.5 km s 1 . 


Table 2. Probable members based on V r . [Fe/H] values marked with * are derived using a V magnitude estimated from 2MASS 
J and K a magnitudes. R = radial distance from cluster centre. Status: m = probable new cluster member star, h = probable 
extratidal cluster halo star. 


2MASS 

ID 

RA 

Dec 

Vr 

(kms -1 ) 

AV r 

(kms -1 ) 

[Fe/H] 

(dex) 

A [Fe/H] 
(dex) 

fL a COs(S) 
(mas yr — - 1 ) 

M<5 

(mas yr — - 1 ) 

R 

(arcmin) 

Status 

05111553-4018102 

5:11:15.54 

-40:18:10.25 

318.7 

2.7 

-1.27 

0.09 

-0.1 

0.7 

36.1 

h 

05135811-3956592 

5:13:58.11 

-39:56:59.25 

324.7 

1.1 

-1.38* 

0.06 

-1.3 

-0.1 

6.0 

m 

05140637-4024514 

5:14:06.37 

-40:24:51.40 

321.3 

1.6 

-1.02* 

0.05 

-8.5 

-2.9 

22.1 

h 

05140736-4000020 

5:14:07.36 

-40:00:02.04 

312.7 

1.6 

-1.15* 

0.05 

-10.8 

13.0 

2.8 

m 

05141174-4002180 

5:14:11.74 

-40:02:18.07 

319.0 

1.6 

-0.90* 

0.03 

- 

- 

1.1 

m 

05141359-4001109 

5:14:13.59 

-40:01:10.98 

314.9 

1.6 

-1.56* 

0.07 

-0.3 

4.7 

2.1 

m 

05141398-4000194 

5:14:13.99 

-40:00:19.39 

321.6 

1.6 

-1.44 

0.06 

-1.7 

12.1 

2.8 

m 

05142006-4005116 

5:14:20.06 

-40:05:11.69 

315.2 

2.2 

-1.19 

0.06 

3.3 

0.3 

3.5 

m 

05142069-3958560 

5:14:20.70 

-39:58:56.02 

311.2 

1.6 

-1.29 

0.05 

4.7 

-0.7 

4.7 

m 

05142088-4005477 

5:14:20.89 

-40:05:47.74 

318.8 

1.3 

-1.21* 

0.05 

-2.1 

2.1 

4.0 

m 

05142994-4000251 

5:14:29.95 

-40:00:25.08 

318.3 

1.3 

-1.08 

0.04 

-4.6 

1.3 

5.0 

m 

05153906-3947485 

5:15:39.06 

-39:47:48.52 

318.6 

1.3 

-1.08 

0.04 

0.2 

-1.9 

23.2 

h 

05162438-4019326 

5:16:24.38 

-40:19:32.62 

317.9 

1.1 

-1.10 

0.04 

-1.5 

-0.4 

31.2 

h 


ing we will refer to stars in the first group as non-members, 
stars in the second group inside the tidal radius as cluster 
stars or members while stars outside the tidal radius will 
be referred to as cluster halo stars. In the second group, we 
recover 106 of 124 of the red giant cluster members in the 
CIO group and the four red giant cluster members observed 
in the YG08 group. We also identify 13 probable new clus¬ 
ter members or cluster halo stars of NGC 1851. The right 
panel shows an expanded plot of the V r distribution, show¬ 
ing all the observed cluster member and cluster halo stars 
and a vertical line showing the V r of NGC 1851. As well 
as the two peaks we also observe a small number of stars 
with V r values intermediate between the two peaks (eight 


between 200 kms -1 and the lower V r membership limit of 
308 kms -1 ) and three stars with very high radial velocities 
of 371.8, 388.8 and 394.8 kms -1 . These stars are discussed 
further in Sect. 8. 

The mean and standard deviation of the V r distribution 
of cluster members and cluster halo stars is 319.3 kms -1 and 
4.4 kms -1 , which is similar to previous estimates for 
NGC 1851 (320.5 ± 0.6 kms -1 from Harris 1996 (2010 
edition). The V r of these stars strongly suggests that all 
these targets are likely cluster members or cluster halo stars, 
and, where available, this is supported by their positions, 
proper motions, metallicities (Sect. 5) and positions on the 
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colour-magnitude diagrams (Sect. 7). Details of these stars 
are listed in Table 2. 

There are no indications of an overdensity at 
~180 kms^ 1 in our observed distribution that might corre¬ 
spond to the tidal tail at this V r reported by Sollima et al. 
(2012). However the number of stars for comparison is small, 
especially if we exclude the number of stars that are pre¬ 
dicted by the Milky Way model generated by the Galaxia 
code distribution to be in this range of V r . Therefore our 
observations do not provide additional support for the exis¬ 
tence of a stellar stream at ~180 kms -1 . 


4 SPATIAL DISTRIBUTION AND PROPER 

MOTIONS 

Figure 4 shows the spatial (RA and Dec) distribution of ob¬ 
served stars. Of the 13 stars classified as probable new clus¬ 
ter members or cluster halo stars on the basis of V r , nine are 
within our adopted Carballo-Bello et al. (2012) tidal radius 
of 11.7 arcmin and are identified as probable current mem¬ 
bers of NGC 1851. A further four are outside at distances up 
to ~3.1 times the adopted tidal radius. Members of a stellar 
halo of NGC 1851 are expected to share its distinctive V r 
signature, so these are identified as probable extratidal clus¬ 
ter halo stars. The identification of these cluster halo stars 
provides further evidence of the existence of a stellar halo 
surrounding NGC 1851. Details of these stars are listed in 
Table 2. 

The absolute proper motions of 106 of the CIO stars, 
nine out of ten of the probable new cluster members and 
all four probable new cluster halo stars of NGC 1851 were 
available in the UCAC4 catalogue (Zacharias et al. 2013). 
The proper motion of NGC 1851 (fi a cos(5 ) = 1.28 ± 
0.68 mas yr _1 , = 2.39 ± 0.65 mas yr -1 ) was obtained 

from Dinescu et al. (1999). Figure 5 shows a plot of the abso¬ 
lute proper motions of the observed CIO stars (only 20 stars 
are plotted to avoid overcrowding) and 12 of the the 13 prob¬ 
able new cluster members or cluster halo stars of NGC 1851. 
The proper motion of NGC 1851 is shown by the intersection 
of the dashed lines. It is interesting to note that 15 of the 
106 CIO stars have UCAC4 catalogue proper motions in RA 
and/or Dec greater than 50 mas yr -1 and 34 have proper 
motions in RA and/or Dec greater than 25 mas yr _1 . This 
scatter in proper motion values may simply reflect the dif¬ 
ficulty of making proper motion measurements in crowded 
stellar fields. The proper motions of the probable new clus¬ 
ter members and cluster halo stars are listed in Table 2 and 
all are consistent with the proper motion of NGC 1851. 


5 METALLICITY FROM THE Ca TRIPLET 
LINES 

Applying the method of Armandroff & Da Costa (1991) we 
used measurements of the equivalent widths ( EW ) of the 
CaT lines to calculate stellar metallicity ([Fe/H]). A number 
of CaT empirical relations have been proposed by different 
authors. They connect a linear combination of the equiva¬ 
lent widths of the CaT lines and the luminosity of the star 
(usually the V magnitude of the star relative to the HB of 
the cluster: V — Vhb ) to its [Fe/H] value. 


We used the IRAF task splot to measure the EW s of 
the CaT lines (8498 A=> EW $498 , 85 42 A=> EWs 54,2 and 
8662 A=> EWs 662 ) in the red channel spectra of the known 
CIO and YG08 members plus the probable new cluster mem¬ 
bers and cluster halo stars. Gaussian profiles were fitted to 
these lines to measure the EWs. 

The calibration procedure is valid for V — Vh b brighter 
than ~0.0 and for a range ~2 in V — Vhb luminosities. 
To include all 13 of the probable new cluster members and 
cluster halo stars, a faint limit of V — Vhb < 0.3 and a 
bright limit of —1.7 < V — Vhb were chosen. These criteria 
were also applied to the CIO and YG08 stars to select a 
comparison sample. 64 of the 106 stars in the CIO group 
matched these criteria, but the YG08 stars are somewhat 
brighter and did not fit within the limits. 

Fig. 6 shows the sum of the measured EWs of the CaT 
lines (-EW 8542 + EWseei) of the 13 new cluster members and 
cluster halo stars plus the 64 CIO stars, plotted against the 
magnitude difference from the horizontal branch V — Vhb- 

To calibrate the EWs in terms of [Fe/H] we adopted the 
empirical calibration equation of Yong et al. (2014): 

[Fe/H] = 0.59(±0.04)W' - 2.81(±0.16) (2) 

where the reduced equivalent width W' is : 

W' = FW8542 + EWsgq2 + 0.60(±0.02)(F - Vhb ) (3) 

V magnitudes were listed in the catalogue for all 64 of 
the CIO stars. V magnitudes were available from APASS 
(Henden et al. 2009) for seven of the 13 probable new clus¬ 
ter members and cluster halo stars. For the remaining six 
stars, an estimate Vj,k„ was made from 2MASS J and K a 
magnitudes. The calibration formula used was that used to 
calculate V magnitudes for the GALAH survey input cata¬ 
logue (De Silva et al. 2015): 

Vj, Kb = K s + 2(J -K a + 0.14) + O.382e ((j-Xa-0 ' 2)/ °' 50) (4) 

There are differences of up to ± ~0.5 mag when V mag¬ 
nitudes obtained with this calibration are compared with lit¬ 
erature values for CIO stars, which translates into an error 
of ± ~0.4 dex in [Fe/H]. Given this, the values of [Fe/H] for 
these stars are best regarded as indicative values only. 

Details of the [Fe/H] of the 13 probable new cluster 
members and cluster halo stars are listed in Table 2, with 
the [Fe/H] of stars obtained using an estimated V magnitude 
marked with an asterisk. 

Fig. 7 shows the [Fe/H] distribution of the 13 probable 
new cluster members and cluster halo stars and the 64 CIO 
members of NGC 1851, with the predicted [Fe/H] distribu¬ 
tion of the Milky Way model generated by the Galaxia code 
(Sharma et al. (2011), see Sect. 6) overplotted as a grey his¬ 
togram. The model distribution is normalised to the number 
of observed known members (64) plus the number of prob¬ 
able new cluster members and cluster halo stars (13). The 
distribution shows a prominent peak with a mean [Fe/H] = 
-1.10 that we identify as members of NGC 1851. This group 
includes the 64 homogeneously selected RGB stars from the 
CIO group and the 13 probable new cluster members or clus¬ 
ter halo stars of NGC 1851. 

The mean and standard deviation of the observed 
[Fe/H] distribution are -1.10 and 0.11 respectively. Seven 
of our probable new cluster members or cluster halo stars 
lie within ±1 <t of the [Fe/H] metallicity of NGC 1851 and 
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Figure 4. Left panel: Spatial distribution of observed stars for the complete 2° field of 2dF/AAOmega. Right panel: Expanded version 
showing the observed stars inside the tidal radius. The green cross and circle show the Harris 1996 (2010 edition) catalogue centre position 
and the tidal radius of 11.7 arcmin. Grey crosses show non members, blue circles denote known CIO members, green squares are known 
YG08 members and red triangles are probable new cluster members or cluster halo stars. 



Figure 5. Absolute proper motions of 20 of the observed CIO 
stars (blue) and 12 of the probable new cluster members or cluster 
halo stars of NGC 1851 (red). The proper motion of NGC 1851 
is shown by the intersection of the dashed lines. 



Figure 6. The sum of the EW s of the CaT lines plotted against 
the magnitude difference from the horizontal branch V — Vhb- 
Blue circles indicate known CIO cluster members, red triangles 
are probable new cluster members and cluster halo stars with V 
magnitudes from APASS and red stars are probable new clus¬ 
ter members and cluster halo stars with V magnitudes estimated 
from J and K magnitudes. 


10 within d=2cr, so most of these stars have metallicities that 
are consistent with that of NGC 1851. 
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Figure 7. [Fe/H] distribution of cluster members and cluster halo 
stars with the predicted distribution of the Milky Way model 
generated by the Galaxia code (Sharma et al. 2011) overplotted 
as a grey histogram. Blue indicates known CIO cluster members 
and red denotes stars that are probable new cluster members and 
cluster halo stars. The vertical dashed line shows the Harris 1996 
(2010 edition) catalogue [Fe/H] of -1.18. 

6 COMPARISON OF OBSERVED V R AND 

[Fe/H] DISTRIBUTIONS WITH A MILKY 

WAY MODEL 

To investigate further whether the observed V T and [Fe/H] 
distributions were associated with NGC 1851, we compared 
our observations with Milky Way models generated by the 
Galaxia code (Sharma et al. 2011). 

The Galaxia code was used to generate a synthetic cat¬ 
alogue of stars covering the same area of sky that can be 
compared with a set of observations. Given colour, mag¬ 
nitude and spatial constraints, Galaxia predicts radial ve¬ 
locities, metallicities and other properties of stars in that 
area of sky. The Besangon Milky Way model (Robin et al. 
2003) was used for the disk and the simulated N-body mod¬ 
els of Bullock & Johnston (2005) for the stellar halo. To 
compare our observations with the model the code was run 
to generate 1000 models covering the 2° field centred on 
the Harris 1996 (2010 edition) catalogue centre position of 
NGC 1851 with no colour or magnitude restrictions. V r and 
[Fe/H] distributions were then constructed after applying 
the same colour and magnitude criteria for selection from 
the 2MASS catalogue that were used in the target selec¬ 
tion process (0.1 < J — H < 0.9, 0 < J — K s < 1.2 and 
13 < J < 14.5). The V r distribution was then summed and 
normalised to have the same number of stars as the number 
of observed non-member stars, while the [Fe/H] distribution 
was normalised to have the same number of stars as the num¬ 
ber of observed known and probable new cluster members 
and cluster halo stars. 

In Fig. 3 and Fig. 7 the observed and model V r and 
[Fe/H] distributions are plotted. It is immediately apparent 


that the cluster stars are not reproduced by the models in ei¬ 
ther the V r or the [Fe/H] distributions. We then performed a 
Kolmogorov-Smirnov (K-S) test for the null hypothesis that 
the model V r and [Fe/H] distributions and the observed V r 
and [Fe/H] distributions are drawn from the same continu¬ 
ous distributions. 

For the V r distribution, this was done in the V r range 
from 250 to 400 kms -1 to compare NGC 1851 member stars 
with the average Milky Way model distribution. The mean 
and standard deviation of the observed V r distribution were 
319.3 kms -1 and 4.5 kms -1 respectively, compared with 
306.7 kms -1 and 40.1 kms -1 for the average model distri¬ 
bution. The p value calculated by comparing the observed 
and model distributions with the K-S test was 3x 10 -35 . For 
[Fe/H] the mean and standard deviation of the observed dis¬ 
tribution were -1.108 and 0.097 respectively, compared with 
-0.309 and 0.373 for the average model distribution. The 
p value calculated by comparing the distributions with the 
K-S test was 9.4x 10 -43 . 

The p values for comparing the observed and model V r 
distributions are effectively zero, hence it is unlikely that 
the observed and predicted V r samples are drawn from the 
same distributions. The Milky Way model generated by the 
Galaxia code predicts just three stars within the same spa¬ 
tial, colour, magnitude and V r ranges as the observed stars. 
Similarly, the p value for the [Fe/H] distribution is effectively 
zero, hence it is unlikely that the observed and predicted 
[Fe/H] samples are drawn from the same distributions. The 
model predicts just six stars within the same spatial, colour, 
magnitude and [Fe/H] ranges as the CIO stars. Most stars 
in the observed distributions are thus either current mem¬ 
bers, probable new cluster members or cluster halo stars of 
NGC 1851. 


7 COLOUR-MAGNITUDE DIAGRAMS 

Fig. 8 shows 2MASS and optical CMDs (J vs J — K s and 
V vs B — V) of target stars. The left panel clearly shows 
the 2MASS selection limits of the third group of target stars 
(0 < J—K s < 1.2 and 13 < J < 14.5). The right panel shows 
the V vs B — V CMD of all the known CIO and YG08 clus¬ 
ter members and probable new cluster members and cluster 
halo stars (with directly measured B and V magnitudes). 
The stars from the photometric study of Walker (1992) are 
also plotted and the split nature of the HB in the cluster is 
obvious. The CIO stars, three of the four YG08 stars and 
seven of the probable new cluster members and cluster halo 
stars that had directly measured B and V magnitudes all 
nicely follow the RGB ridge line of the cluster population. 


8 HIGH V R STARS 

There are eight stars with V r between 200 kms -1 and 
the V r of NGC 1851 as listed in Table 3. These stars are 
likely not cluster members so their distance and luminos¬ 
ity are unknown. Any V — Vhb measurements would be 
invalid, so metallicities could not be calculated for these 
stars. Proper motions obtained from the UCAC4 catalogue 
(Zacharias et al. 2013) are also tabulated. The proper mo¬ 
tion of NGC 1851 is p a cos(5) = 1.28 ± 0.68 mas yr _1 , 
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Figure 8. Left panel: 2MASS (J vs J — K s ) CMD. The grey crosses show observed non member stars, blue circles indicate known CIO 
members, green squares known YG08 members and red triangles denote stars that are probable new cluster members and cluster halo 
stars. Right panel: V vs B — V CMD. Here the grey crosses show Walker (1992) stars, other symbols are the same as for the left panel 
except that only the seven probable new cluster members and cluster halo stars with directly measured B and V magnitudes are plotted. 


Table 3. Other stars with high V r . Stars above the line are those with V r between 200 kms" 1 and the lower V T membership 
limit for NGC 1851 of 308 kms -1 , those below the line have V r greater than the V r of NGC 1851. R. = radial distance 
from cluster centre. 


2MASS 

RA 

Dec 

V r 

AV r 

fioiCos(8) 

ns 

R 

identification 



(kms -1 ) 

(kms -1 ) 

(mas yr — - 1 ) 

(mas yr — - 1 ) 

(arcmin) 

05122418-4051595 

5:12:24.18 

-40:51:59.53 

210.3 

1.6 

9.3 

-18.4 

52.9 

05130199-3931439 

5:13:02.00 

-39:31:43.90 

283.6 

1.2 

13.7 

-1.8 

33.5 

05134950-4013007 

5:13:49.50 

-40:13:00.77 

201.5 

1.1 

35.2 

-23.8 

10.7 

05141567-4001444 

5:14:15.68 

-40:01:44.39 

301.1 

1.6 

5.9 

6.2 

2.0 

05152527-4012529 

5:15:25.27 

-40:12:52.94 

243.7 

1.3 

16.4 

8.7 

18.1 

05162712-3931579 

5:16:27.12 

-39:31:57.95 

213.6 

0.9 

-1.9 

8.3 

41.0 

05163977-4039434 

5:16:39.77 

-40:39:43.39 

288.0 

2.2 

60.4 

-0.9 

47.0 

05182129-4008551 

5:18:21.30 

-40:08:55.15 

292.4 

1.3 

1.8 

3.7 

49.1 

05133314-4100160 

5:13:33.14 

-41:00:16.05 

371.8 

1.7 

-1.0 

5.4 

57.8 

05135111-3958381 

5:13:51.11 

-39:58:38.13 

388.8 

1.3 

5.4 

-2.0 

5.1 

05184308-4027077 

5:18:43.09 

-40:27:07.78 

394.8 

1.5 

4.2 

3.4 

58.1 


/is = 2.39 ± 0.65 mas yr^ 1 ) (Dinescu et al. 1999). Sev¬ 
eral stars have proper motions that are quite different, but 
as mentioned in Sect 4 there are also significant numbers of 
CIO stars in the UCAC4 catalogue with proper motions that 
are considerably different from that of NGC 1851. One star, 
2MASS 05141567-4001444, is located just 2 arcmin from the 
cluster centre, has a proper motion similar to NGC 1851 and 
a V r of 301.1 kms' 1 which is just outside the V r criteria we 
adopted for membership, so it could conceivably be a clus¬ 
ter member. If cluster membership is assumed, the [Fe/H] 
calculated using the methods of Sect 5 using an estimated 
Vj,k s from 2MASS J and K s magnitudes, is -0.41. This is 
high for a cluster member, even allowing for errors in calcu¬ 
lating [Fe/H] with estimated V magnitudes (discussed above 
in Sect 5). Therefore this star was not included in the list of 
probable members in Table 2. 

There are also three stars with V r ~50 to 


75 kms 1 higher than that of NGC 1851 (05133314-4100160 
at 371.8, 05135111-3958381 at 388.8 and 05184308-4027077 
at 394.8 kms' 1 ) listed in Table 3. As these stars are 
also likely not cluster members, metallicities could not 
be calculated for them. Based on the similarity of their 
Vr and proper motions we speculate that the stars might 
be associated. None of the GCs that are also possibly 
associated with the Canis Major dwarf galaxy (NGC 1904, 
NGC 2298 and NGC 2808) have a V r that is consistent with 
these stars, in fact there are no GCs in that area of the 
sky that have a similar V r . Without further information, at 
this stage we tentatively identify all three as high velocity 
Galactic halo stars. 
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9 SUMMARY AND CONCLUSIONS 

We report the first results of a spectroscopic survey of 
NGC 1851 and its surroundings, with good quality spec¬ 
tra of 1149 unique stars. Based on V r obtained from Fourier 
cross-correlation of the wavelength range around the near 
infrared calcium triplet absorption lines, we identify 

(1) nine probable new cluster members inside the tidal 
radius 

(2) four probable extratidal cluster halo stars at distances 
up to ~3.1 times the tidal radius 

(3) 106 red giant cluster stars in the Carretta et al. (2010, 
2011 ) study 

(4) four red giant cluster stars observed in the 
Yong & Grundahl (2008) study 

(5) eight stars with V T values intermediate between those 
of the Galactic held stars and NGC 1851 

(6) three stars with very high radial velocities of 371.8, 
388.8 and 394.8 kms -1 , and 

(7) 1015 non-member Galactic disc or halo held stars with 
-100 kms -1 < V r < 200 kms -1 

Proper motions were available for eight of the nine prob¬ 
able new cluster members and the four probable extrati¬ 
dal cluster halo stars and all are consistent with that of 
NGC 1851. The four extratidal cluster halo stars provide fur¬ 
ther conhrmation of the hndings of Olszewski et al. (2009), 
Carballo-Bello & Martfnez-Delgado (2010), Marino et al. 
(2014) and Kunder et al. (2014) that NGC 1851 is sur¬ 
rounded by a stellar halo system and does not have a clas¬ 
sical tidally-limited prohle. Apart from the cluster and halo 
stars our observed radial velocity distribution agrees with 
the expected distribution from the Besangon disk/N-body 
stellar halo Milky Way model generated by the Galaxia code, 
and the probability that the V r distribution of the cluster 
and halo stars are drawn from the same distribution as the 
model is effectively zero. In particular there are no indi¬ 
cations of an overdensity at ~180 kms -1 in our observed 
distribution that might correspond to the tidal tail at this 
V r reported by Sollima et al. (2012), however the number of 
observed stars in our sample with comparable V r is small. 

The metallicities of the Carretta et al. (2010, 2011) 
stars and the probable new cluster members and cluster halo 
stars were estimated using equivalent width measurements 
of the calcium triplet absorption lines. Within the limita¬ 
tions of this method discussed above, the metallicities of 
most of the probable new cluster members and cluster halo 
stars were found to be consistent with that of NGC 1851. 
As with the V r distribution, the probability that the [Fe/H] 
distribution of the cluster and halo stars are drawn from the 
same distribution as the model is effectively zero. 

We also identified three stars with unusually high V T , 
~50 to 75 kms -1 higher than that of NGC 1851. As the 
velocity separation is large it is likely that they are not as¬ 
sociated with NGC 1851. Based on the similarity of their 
V r and proper motions we speculate that the stars might 
be associated and tentatively identify them as high velocity 
Galactic halo stars. 

In our second paper (in preparation) we will address the 
determination of C and N abundances using the blue channel 
of the spectra and discuss the questions of the relationship 
of the C-N and O-Na anticorrelations, whether there is a 


scatter in the sum of the C+N+O abundances and whether 
the sum correlates with metallicities or heavy-element abun¬ 
dances. 
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